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It is shown that the quasi-localized states in weakly dis- 
ordered systems can lead to the non-analytical distribution 
of level curvatures. In 2D systems the distribution function 
P(K) has a branching point at K = 0. In quasi-lD systems 
the non-analyticity at K — is very weak, and in 3D metals 
it is absent at all. Such a behavior confirms the conjecture 
that the branching at K = is due to the multi-fractality 
of wave functions and thus is a generic feature of all critical 
eigenstates. The relationsip between the branching power and 
the multi-fractality exponent r/(2) is derived. 
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The problem of energy level statistics in complex quan- 
tum systems has recently attracted much attention. It is 
a key problem for the description of mesoscopioc fluc- 
tuations in small but macroscopic disordered electronic 
systems at low temperatures 1 . It is also of great fun- 
damental importance, since level statistics is a proper 
language to describe the signature of chaos in a quantum 
system whose classical counterpart is chaotic. 

The most of investigations in the field has been aimed 
to describe the spectral statistics themselves indepen- 
dently of the nature of wave functions. Nonetheless, it 
is the space structure of wave functions that determines 
the properties of spectral statistics. In systems where 
the wave functions are extended and essentially struc- 
tureless (e.g. in good metals) the predictions of the ran- 
dom matrix theory (RMT) 2 is proved to work well 3 . 
The corresponding Wigner-Dyson (WD) spectral statis- 
tics is characterized by the level repulsion. On the other 
hand, in the region of localized states, energy levels are 
almost uncorrelated obeying the Poisson statistics. The 
spectral statistics at the mobility edge 4,5 proves to be 
universal but different from both the Poisson and the 
WD statistics. Since the wave functions at the mobility 
edge are still extended, the energy levels repel each other 
at small energy separations u> of the order of the mean 
level spacing A, the correlations decaying as a power law 
5 for u; ^> A. However, the typical fluctuation of the 
number of levels SN in an energy window of the width 
E = NA 3> A is much bigger than the RMT prediction 
and behaves in a quasi-Poisson way 6 ' 7 : 



((SN) 2 ) = X N^>\nN. 



(1) 



tions. A typical wave function at the mobility edge can 
be considered as a superposition of peaks with the broad 
distribution of heights and widths. The corresponding 
quantum- mechanical probability density |^(r)| 2 is sparse 
and is not negligible only in a small fraction of volume 
that scales with the system size L as I D , D = d — rj 
being the fractal dimensionality in the d-dimensional 
space. However, two wave functions of different ener- 
gies \E — E'\ A overlap well 8 with each other since 
they are strongly correlated in a given point: 



\* E (r)\P\* E ,(r)\P) oc \E - E'\- v(p)/d - 
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Such slowly decaying correlations are typical of the crit- 
ical wave functions. Both in good metal and in the in- 
sulator wave functions of different energies are locally 
uncorrelated. The statistics of wave functions expressed 
by Eq.(^|) exhibits multi-fractality 9 , since the fractal di- 
mensionality D(p) = d — 7]{p) depends on the power p. 
The main effect of multi-fractality is described by the ex- 
ponent 77(2) = T). In particular the level compressibility 
X in Eq.(|l|) has been found 7 to depend on 77 as follows: 
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Such a complex behavior of the critical spectral statistics 
is due to the complex statistics of the critical wave func- 



Multi-fractality of wave functions is typical not only for 
critical states at the mobility edge. Weak multi-fractality 
holds 10-12 for two-dimensional disordered systems in the 
regime of weak localization where the system size L is 
much less than the localization radius Rq. 

In general, multi-fractality can be considered as a re- 
sult of proliferation of the quasi-localized states. The lat- 
ter notion 13,14 means rare states in disordered systems 
where ^(r)! 2 has a sharp peak on top of the structureless 
extended background (in 3D metals) or localized states 
with anomalously small localization radius (for d < 2). 
The existence of such states exhibits itself in the slow 
relaxation of current first discovered in Ref. [11] using 
the complicated renormalization-group analysis of an ex- 
tended nonlinear sigma-model in 2 + e dimensions. 

Very recently the problem of current relaxation in dis- 
ordered conductors has been reconsidered 14 ' 15 by an el- 
egant instanton approximation 14 applied to supersym- 
metric version of the nonlinear sigma-model 3 . In these 
works the main result of the previous consideration 11 
has been confirmed for 2D systems. However, the new 
method was able to describe some unknown regimes of 
current relaxation and to set correct limits of validity 
for the regimes found earlier. Later the same idea 14 
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has been applied 12 to find directly the distribution of 
^(r)! 2 and the distribution of local densities of states 
16 . Thus the existence of the quasi-localized states has 
been proved and the corresponding configuration of the 
random impurity potential has been found 14 ' 17 . 

The question we address in this Letter is whether and 
how the quasi-localized states and multi-fractality man- 
ifest itself in level statistics other than the level number 
variance. The natural candidate is the statistics of level 
curvatures K n first introduced by Thouless 18 : 



1 d 2 E n {4>) 
n ~ A d<b 2 



(4) 



Here E n (<p) is an energy level depending on the twist 
of phase in the boundary condition ty(r±,x = 0) = 
^(r^, x — L) e 1 ^ . Such a phase shift arises if an external 
magnetic flux threads a sample of the toroidal form. 

The main idea of Ref. [18] is that it is possible to distin- 
guish between the localized and extended states studying 
the sensitivity of spectrum to the twist of phase in the 
boundary conditions. This sensitivity is maximal for ex- 
tended states and negligible for the localized ones. It 
is clear that the existence of the quasi-localized states 
should lead to an enhancement in the probability to find 
a state with a small value of the random quantity K n . 
Thus the quantity of interest is the deviation of the dis- 
tribution function P(K) from the form 19 ' 20 that follows 
from the random matrix theory. 

This deviation has been recently studied numerically 
21 and analytically 20,22 using an expansion in powers of 
the inverse dimensionless conductance g . The main 
result of the numerical simulations 21 on the mobility 
edge of the 3D Anderson model is that the deviation 
5P(K) = P{K) - Pwd(K) cx \K\ 2 - a has a branching 
point at | .fiT | =0 with a 0.4. Although the analytical 
result 22 describes correctly the sign and shape of the 
deviation SP(K), it is unable to explain the branching at 
K = which should be a non-perturbative in 1/g effect. 
It has been also conjectured in 21 that the branching is 
due to the multi-fractality of critical wave functions and 
thus the new critical exponent a should be related with 
the multi-fractality exponent 77. 

In this Letter we will show using the instanton approx- 
imation similar to that used in 14 that the level curva- 
ture distribution P{K) has, indeed, a branching point at 
if = in 2D systems under the weak-localization con- 
dition i?o 3> L. In 3D metals and quasi-lD systems 
with Rq L the distribution function deviates from the 
WD form but all its derivatives are finite at K = 0. 
This result favors the statement 21 that the branching 
in P(K) at K = and multi-fractality are related with 
each other, since for good conductors the multi-fractality 
is also present only for d=2. 

We have found an explicit relationship between a and 
r) that we believe is valid at the mobility edge near the 
Anderson transition too: 



Another relationship can be found using Eq.(||): 

a = 3-~x~ 1 , (2/d > X > 2/3d). (6) 

The first of these equations relates the parametric level 
statistics P(K) with the multi-fractal statistics, Eq.(||), 
of critical wave functions. The second one establishes a 
connection between the parametric and the conventional 
level statistics at the mobility edge. 

The level curvature distribution function P{K) = 
(S(K - K n )5{E - E n )) can be expressed 20 ' 22 in terms of 
the supersymmetric nonlinear sigma-model derived in 3 
from the model of non-interacting electrons in a random 
Gaussian potential: 



P(K) = lim ( 

0^0 



3? / VQ A[Q] e 



-F[Q;K,<f>] 



(7) 



a = 3-4/77, (77 > 4/3). 



(5) 



where A[Q] is a certain pre-exponent that is irrelevant 
for the present consideration and 

F[Q- K, 4>] = ^g J Str {dQ - i<f>n [Q, f]) 2 d d r + (8) 

+ i^K<j? J Str[kQ] d d r. 

Here g = A~ 1 D/L 2 is the dimensionless conductance, D 
is the diffusion coefficient; the unit vector n is directed 
along the x-axis. The field Q(r) is a 8 x 8 super-matrix of 
the certain symmetry 3 that obeys the condition Q 2 = 1 . 
The symbol Str[] stands for the super-trace. The matri- 
ces t = er z (l + A)/2 and A = diag(l; —1), describe break- 
ing of time-reversal symmetry and the symmetry between 
retarded and advanced Green's functions, respectively. 
Here and below we use the dimensionless coordinates- 
ordinates r = {r^}, where —1/2 < n < 1/2. 

One can see that the level curvature K enters Eq.(||) 
in the same way as the frequency uj in the problem of 
long-time current relaxation 14 . Therefore from the re- 
sults of Ref. [14] one would expect the instanton approx- 
imation to work well for the Fourier-transform P(X) = 
J dK exp[—iKX] for A ^> 1 rather than for the distribu- 
tion function P{K) itself. 

It is very important that in order to find the distri- 
bution function P(K) or P(X) one should do the limit 
<\> — > 0. In this limit a nonzero contribution to the func- 
tional F[Q; K, <f> — ► 0] can arise only from the components 
in the non-compact boson-boson sector of the field Q that 
can take arbitrary large values. 

Another simplification comes from the fact that inte- 
gration over Grassman variables requires an expansion 
of exp{— F[Q; K, </>]} in powers of the Grassman compo- 
nents of the field Q. Therefore such components can 
only contribute to the pre-exponent A in the expression 
P(X) = Ae~ s( - X \ where the effective action 5(A) > 1. 
The main, exponential effect is due to the conventional 
variables that parametrize the matrix field Q. Thus using 
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the Efetov's parametrization 3 for the orthogonal ensem- 
ble we have Q — V~ 4 HV ', where H contains the "non- 
compact" angles #1,2(1") > 1 and V contains "phase" vari- 
ables (f(r) £ [0, 2tt] and x( r ) G [0, 2ir]. All functions must 
obey the periodic boundary conditions. 

By varying the functional F[Q;K,<p] + iKX over these 
functions we find the saddle-point equations to be sat- 
isfied by the choice d\ = 0, 9\ = 62 = 9 and <p(r) = 
4>v(t) ^ 0. The equations for the rest two functions read: 



d 2 9 + <j) 2 [K - (dv - n) 2 ] sinhtf = 0, 



d[(dv-n) (cosh 9-1)] = 0. 



(10) 



The variation over K = —2ign leads to the self- 
consistency equation, Eq.(p^). 

The limit — > is done simply by absorbing (f> 2 into 
9. We introduce 6 = 9 + ln0 2 . Then in the limit 0^0 



we have sinh0 = cosh# = |e <f> 
from Eqs.(|),©. 

One can solve Eq.(|lC|) to obtain: 

(dv — n) = rot A. e~' 



3> 1 so that 4> drops 



(11) 



where A is an arbitrary vector function. In what follows 
we limit ourselves by considering rotA. = —n/N = const. 
Integrating Eq.(p"l|) over space and using the periodic 
boundary conditions for v(r ) we obtain the second self- 
consistency equation Eq.(O). 

All equations take a more symmetric form if we make 
a shift 9 = u- C, where sinhC = JVk" 1/2 (k - iV~ 2 )/2: 



d 2 



7 sinhw = 0, 



A= 7 



A = 



872 AT2 



d d 7 



e u cfV, 



(12) 



(13) 



(14) 



where 7 2 = y/R/N. 

Solving these equations with the periodic boundary 
conditions one finds u(r, A), N(X) and 7(A) that enter 
the instanton action 5(A) = F[Q; K, 0-^0]+ iKX\ s : 



- lnP(A) ~ 5(A) = -g I (du) 2 d d p + 2gj N X. (15) 

Eqs . (p~2|) , (|T^) are analogous to those derived in 14,15 for 
the problem of current relaxation in the orthogonal en- 
semble. The second self-consistency equation (|13| ) is typi- 
cal for our problem which is essentially a crossover prob- 
lem from the orthogonal to the unitary ensemble. An- 
other important difference is that we are looking for the 
periodic solutions to the above equations with the period 
a = 1 rather than the solutions 14 that vanish at the 
boundary of the sample. 



Since 7 2 cx K x l 2 we expect the saddle point to corre- 
spond to 7 2 <C 1 at A — > 00. In this case the nonlin- 
ear term in Eq. (TL2) can be approximated by 7 2 sinh u « 

2 I I — 
-fe-el"l sign(u). Thus we come to the Liouville equation. 

Let us first consider the case of a quasi- ID sample. 
Then the proper periodic solution to the Liouville equa- 
tion takes the form: 



4fc 2 7 2 cosh 2 (kx) 



(\x\ < 1/4), (16) 



(9) where 4fc 2 = 7 2 cosh 2 (fc/4). In order to obtain the so- 
lution u(x) for 1/4 < \x\ < 1/2 one should simply con- 
tinue the solution, Eq.jjlq), anti-symmetrically about the 
points x — ±1/4. Because of such anti-symmetricity we 
have 



d d r = / e +u d d r 



(17) 



Finally for the case of quasi-lD sample we obtain N 1 
8k « 81n(l/ 7 4 ), 7 4 = 64nk 3 /X, and 



P(X) = Aexp 



-^ln 2 A 



(18) 



where g\ = 2irg. This result is exactly of the same form as 
the long tails of current relaxation 14 (in the orthogonal 
ensemble) if one replaces t — * A. The region of validity 
1«A< e L l l of Eq.(18) is found from the condition of 
the breakdown of the diffusive regime 14,15 \du(v)\ max = 
L/l, where I is the elastic scattering length. 

In the 2D case we will look for the continuous periodic 
solution to the Liouville equation with the symmetry of 
a square. It is positive inside the square Q with vertices 
at (0; ±1/2) and (±1/2; 0) and is negative in the rest half 
of the sample, being anti-symmetric with respect to all 
sides of the square and thus obeying Eq.(17). 



The problem of finding such a solution to the Liouville 
equation can be solved exactly using the conformal trans- 
formations of a unit circle onto the square f2. However, 
for our purposes it is sufficient to know the solution for 
r <C 1 where it depends only on r. The solution that 
leads to the finite action 5(A) is given by: 



3< r ) = 16b(~/ 2 + br 2 y 



(19) 



where b = Air 4 . Using this solution we find N 1 = I671", 
4 = 8 3 tt 4 /A, and 
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P(X) — AX~ 292 . 



(20) 



where 52 = 2ir 2 g. Besides the validity condition 1 <C 
A <C (L/l) 4 , this result has again exactly the same func- 
tional form as the long tails of current relaxation 14 in 
2D systems. 

In the 3D metal phase, the saddle-point solution for 
P(X) that depends only on one parameter, the dimen- 
sionless conductance, is absent. The reason is that the 
3D Liouville equation does not have a solution that cor- 
responds to a finite action 5(A) in the limit L/l — > 00. 
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It is remarkable that only in the 2D case the instan- 
ton contribution to the Fourier-transform P(A) of P(K) 
has a power-law asymptotic behavior at A ^> 1. In the 
limit 23 L/l — > oo this behavior is not restricted from 
above and the distribution function P{K) must have a 

branching point at A" = 0. Namely, there should exist 

i, - 

a set of critical values of g 2 



(n) 



for 



- < g 2 < 
at \K\ -> 



gn+l 



such as 



the derivative d 2n P{K)/dK 2 



0. The value of a„ is given by: 



a„ = (2n+l)-2ga. 



(21) 



In particular, for g 2 < i, the distribution function P{K) 
should be divergent at K — * 0. For | < g 2 < § it 
should behave like P(K) - P(0) cx |A \ 2 ~ a where a = 
«! = 3 — 2g 2 . One can eliminate the dependence on g 2 
using the expression 12 for the multi-fractality exponent 
r\ = ~. Then we arrive at the result, Eq.(|5|). 

In the quasi-lD case the function P(X) falls off fast 
enough, and all derivatives of P(K) are finite at K = 0: 



P 2n = d 2n P/dK 



2n I 



I K=0 



exp[(2n+l) 2 /2 ffl 



(22) 



In this case we have much softer non-analyticity at K = 
than that in 2D case. However, the function P(K) 
is still not regular at K = 0, since the Taylor series 
p2n A 2 ™/(2n)! has zero radius of convergence. 

Only in the 3D case the quasi-localized states are not 
that important in the metal phase and the correction 
P{K) — Pwd(K) can be found perturbatively in 1/g as 
it has been done in Ref. [22]. 

The above treatment does not allow to take into ac- 
count localization effects explicitly. Thus we are not able 
to derive the level curvature distribution at the mobility 
edge of 3D disordered systems. However, we believe that 
the distribution function at the mobility edge is similar to 
the one in 2D systems, since both cases are characterized 
by the multi-fractal statistics of wave functions. That is 
why we believe the relationship, Eq.(^), that we derived 
for 2D systems in the weak localization regime, to hold 
also in the critical region near the Anderson transition. 

In summary, we have shown that the quasi-localized 
states and multi-fractality of wave functions can lead to 
a non-analytical distribution of level curvatures K. The 
distribution function P(K) is shown to have a branch- 
ing point at K = for 2D systems and a softer non- 
analyticity in quasi-lD systems in the weak-localization 
regime. In the 3D metal the distribution function is 
analytical at K = 0. We argue that the branching in 
P(K) is typical for all critical states where the statis- 
tics of eigenf unctions exhibit multi-fractality. This al- 
lows to conjecture a general relationship between the 
branching power and the multi-fractality exponent rj. 
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The limit L/l — > oo should be taken at a constant g 
by a simultaneous increase in the number of transverse 
channels. 



4 



